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ABSTRACT 


Current understanding of a double vortex thunderstorm involves 
a pair of contra-rotating vortices that exists in the dynamic updraft. 
The pair is believed to be a result of a blocking effect which may 
occur when a cylindrical thermal updraft of a thunderstorm protrudes 
into the upper level air and there is a large amount of vertical wind 
shear between the low level and upper level air layers. 

Eagleman developed a numerical tornado prediction scheme based 
on the double vortex thunderstorm. The Energy-Shear Index (ESI) is 
part of the scheme and is calculated from radiosonde measurements. The 
ESI incorporates parameters representative of thermal instability and 
blocking effect, and indicates appropriate environments for which the 
development of double vortex thunderstorms is likely. 

The ESI and modifications of it were tested using data derived 
from NASA's Fourth Atmospheric Variability Experiment (AVE IV). The 
results showed that the index has some general usefulness in fore- 
casting tornadic outbreaks over a large area, but is probably not 
definitive enough for operational use. At times the results indicated 
large areas of expectant tornadic activity where little actually 
occurred. This deficiency led the author to investigate the 
possibility of incorporating into the index an additional dimensionless 
parameter to take into account the contribution to tornadic development 
associated with the fluid mechanics and shear instabilities of the 
double vortex thunderstorm itself. 


To demonstrate the feasibility of obtaining information on this 
parameter, a wind tunnel investigation was conducted. The cylinder and 
jet in a laminar crossflow at conditions represented by a Reynolds 
number (based on cylinder diameter, crossflow velocity, and kinematic 
viscosity) near 5,000 were investigated. The method of investigation 
was flow visualization. Transition waves in the region extending from 
the cylinder to near three radii downstream of the cylinder's center 
were photographed. 

Since complete simulation of the atmosphere with this wind 
tunnel facility is not possible, the results of the investigation 
cannot be directly related to a physical description of a double vortex 
thunderstorm. The direction for future investigation, however, is 
suggested from this study. 
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CHAPTER I 


NUMERIC AND FLUID DYNAMIC REPRESENTATION OF 
TORNADIC DOUBLE VORTEX THUNDERSTORMS 

There are several methods by which to approach the problem of 
tornado prediction, One method is the use of a numerical prediction 
index, such as the one developed by Eagleman [1]J 

Eagleman suggests that conditions necessary for the production 
of tornadic vortices depend on the development of a double vortex 
thunderstorm. Figure 1.1 is a schematic of the double vortex thunder- 
storm showing the relationship between the thermal updraft and the 
dynamic updraft, Current understanding of a double vortex thunderstorm 
involves a pair of contra-rotating vortices In the dynamic updraft 
which is believed to be the result of a blocking effect between the 
thermal updraft and the upper level air layers. A pair of contra- 
rotating vortices in a tornadic thunderstorm has been fairly well 
substantiated by dual -Doppler radar measurements [2]. 

The Eagleman model of a double vortex thunderstorm is the basis 
for a tornado prediction scheme that incorporates thermal instability 
(the Energy Index, El) and a blocking effect (the Shear Index, SI), 
These two indices are combined when calculating the Energy-Shear Index 
(ESI). The Energy-Shear Index may be computed from data obtained from 


^Numbers in brackets refer to similarly numbered references in 
the Bibliography. 
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INFLOW 


Figure 1.1. Schematic of the double vortex thunderstorm [5] 
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radiosonde measurements, it Is Indicative of appropriate environments 
for development of double vortex thunderstorms but It Is not foolproof. 

The present Investigator tested Eagleman's model, and modifi- 
cations of It, using data derived from NASA's Fourth Atmospheric 
Variability Experiment (AVE IV). NASA conducted the experiment on 
April 24-25, 1975 [3], in which 42 radiosonde stations participated. 

I. THE ENERGY INDEX 

The Energy Index, as developed by Darkow [4], reflects atmos- 
pheric instability. A measure of the instability of the atmosphere is 
the difference of total specific energy between the lower and upper 
level layers. The total specific energy of a unit air mass is the sum 
of the specific enthalpy, gravitational potential energy, latent heat, 
and kinetic energy, 

Ej ■ (c p T) + (g Z) + (L q) + (V 2 / 2) . (1,1) 

Darkow explains that the kinetic energy term is usually two 
orders of magnitude smaller than the other terms. Due to thR 
inaccuracy in reporting upper level humidity, the specific humidity (q) 
is found to be approximately equal to the mixing ratio (W). The latent 
heat of vaporization is assumed to be the value of 0°C. 

Darkow compared average total specific energy profiles of 
environments near tornado touchdowns. The average total specific 
energy profile of 27 proximity soundings taken within 50 miles and 
1 hour and 45 minutes of a tornado touchdown was obtained. Figure 1.2 

shows the profile. The minimum total specific energy of soundings 
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Figure 1.2. Average total specific energy profiles of 
27 tornado proximity soundings [5], 


1 1 1 r 


/ 

\ 

\ 


/ 


/ 


/ 


/ 


/ 


\ Tornado Proximity 
Soundings 

\ 

\ 



cvmn 

OF POOL* 


t , , • * * 


■« * 



taken In proximity of tornadoes occurs near the 500-mb pressure level. 
This fact was used as the Oasis for Darkow's Energy Index (El). 

The El is defined as the difference in total specific energy 
between the 500- and 850-mb pressure level, The total specific energy 
o f a pressure level is calculated by: 


* [t( 0 K) x C (J/gm °K) 

^ r 


g(m/sec 2 ) x z(m) 


L 0 (J/gm) x W(gin/kg) y^ 0 k 0 9 


gm _ 


ihZ) 


E t = T + 9.8 z X ICf 3 + 2.5 W . 


(1.3) 


The Energy Index reflects instability, and when expressed in 
units of calorie per gram (cal/gm), is computed by: 


E t — E t 
1 500 1 850 

“ 4.18 (J/cal } 


(1.4) 


The more negative the El, the more the air mass is expeptpri to 
be unstable. This is because low level air has more energy than upper 
level air. 


II. THE SHEAR INDEX 

The Shear Index (SI) was developed as an indicator of 
appropriate environment for the development of double vortex thunder- 
storms. Severe thunderstorms commonly travel considerably slower than 
the mean environmental wind. The intrusion of warm, moist air into 
upper level air and large amounts of vertical wind shear results in a 
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blocking effect. It was this concept of a blocking effect which led 
Eagleman to the development of a theoretical double vortex model of a 
tornadic thunderstorm. 

The Shear Index reflects changes in wind velocities between low 
level and upper level layers. For a given storm velocity, the SI is 
the number of consecutive upper level layers, I50-mb- thick and taken 
every 50 mb, that have a component of velocity approximately equal to 
and opposing the average velocity of the low level wind. Figure 1.3 
shows the relationship of the six upper level layers to the low level 
layer of the model. The average velocity of an upper level layer is 
computed using velocities measured 25-mb apart throughout the 150-mb- 
thick layer. The average velocity of the low level layer is computed 
using velocities taken 25-mb apart from the surface to the 850-mb level. 
Figure 1.4 shows the vector relationships of relative storm velocities 
and wind velocities. 

Thunderstorm velocities have been observed to vary from 50 to 
80 percent in magnitude and 60 deg in direction either side of the 
average cloud layer wind. Figure 1.5 shows 182 trail storm velocities 
60 deg either side and 50 to 80 percent of the average cloud layer 
velocity. The trail storm velocity of 30 deg left and 55 percent of 
the average cloud layer wind is indicated on Figure 1.5. Shear Indexes 
are calculated for each of the 182 trail storm velocities. The 
velocity for which the maximum SI is calculated is the presumed real 
storm velocity. The average cloud layer velocity used in the calcu- 
lation is the average of the velocities at the 850-, 700-, 500-, and 
300-mb pressure levels. The greater the SI, the greater is the 
expected potential for double vortex formation in thunderstorms. 
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Figure 1.4. Vector relations of storm velocity and wind velocities [5]. 
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III. THE ENERGY-SHEAR INDEX 

To account for thermal Instability and environmental winds the 
Energy-Shear Index (ESI) incorporates the SI and El. The ESI is the 
best empirical linear combination of the SI and El based on a study of 
59 soundings associated with tornadoes [1], A proximity sounding is a 
sounding in the warm air ahead of a cloud front less than 120 miles 
from and within two hours before or one-half hour after a confirmed 
tornado touchdown. A non-proximity sounding is for the same time 
period but more than 200 miles away from a tornado. A precedence 
sounding is one that is greater than 200 miles away from a tornado and 
less than two hours before or one-half hour after a confirmed tornado 
touchdown, Figure 1.(5 displays the 59 soundings and the line of best 
separation between proximity and non-proximity soundings. 

The equation for this line is: 

El = 1/2 SI - 2 (1.5) 

or 

2 El - SI + 4 = 0 . (1.6) 

The ESI is calculated by: 

ESI * 4 - SI + 2 El . (1.7) 

When the ESI is negative, it represents conditions below the line of 
separation and double vortex thunderstorms are expected. For a more 
detailed explanation of the ESI, see Reference [1]. 
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Forecasting Abilit y 

The forecasting ability of the Energy -Shear Index (see 
Appendix A) was tested by Ey [5] using data obtained during Nasa's 
Fourth Atmospheric Variability Experiment (AVE IV). The AVE IV experi- 
ment occurred 0:00 GMT April 24 to 12:00 GMT April 25, 1975. Forty-two 
weather stations in central and eastern United States participated. 

The synoptic conditions at the time of the experiment (by Ey [5]) are 
given in Appendix B. The severe events occurring during the experiment 
(by Turner [6]) are listed in Appendix C. 

The results of Ey's test were deficient for operational use 
because large areas were predicted to have tornadic outbreak, yet 
little activity occurred. In an attempt to reduce the large area that 
falsely predicts tornadic outbreak ("false alarm" area) the present 
investigator modified and tested Eagleman's program (Program 1) three 
times , 

Each modification of Eagleman's program was made to reduce the 
predicted "false alarm" area and expose layers which most affect the 
forecasting ability. The program was modified by altering the pressure 
levels used in obtaining the average wind velocity of the low level and 
the upper level air layers. Tests of the modifications were executed 
on The University of Tennessee IBM 360 computer by a program described 
in Reference [7], 

The surface layer air is well mixed, resulting from interaction 
with the ground and extends up to the cloud base. There, it was con- 
cluded that a better representation of the low level layer velocity 
could be obtained by the average of velocities from the surface to the 
lifting condensation level. 
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The first modification (Program ?.) obtains the low level wind 
from the average of the winds from the surface to the lifting conden- 
sation level, in 25-mb increments. For a dry adiabatic lapse rate of 
10°K/km, the lifting condensation level was calculated by: 


LCL(m) - [T SFS - T D p(°K) 


x 


1 (km) 1,000 (m) 

miTw 



(m) . 


(1.8) 


The predicting ability of Program 2 was tested using AVE IV 
data and compared with the predicting ability of Eagleman's program 
(Program 1), The results of the test are shown in Figures 1.7, 1,8, 
and 1.9 for time periods 24 April 1800 GMT, 24 April 2400 GMT, and 
26 April 0600 GMT, respectively. The -1.0 isopleth is the expected 
line of separation between tornadic and non-tornadic storms. This 
isopleth is shown on the figures for each program. The circled number 
of the program indicates the isopleth for that program. Generally, 
altering the program in this manner did not change the predicting 
ability. Station 213, located near the southeastern coast, for the 
time period 24 April 2400 GMT, was added to the number of stations 
expecting tornadic outbreak where none occurred. The modification of 
the program involving the low level air layer had little effect on 
forecasting. 

A second modification was made to the program to reduce the 
false alarm area. The modification this time involved the number and 
location of upper air layers used when determining blocking potential. 
The number of upper air layers was increased from 6 to 12 layers. The 
location of the additional layers was chosen so as to reflect blocking 
at greater altitudes. The 150-mb-thick layers were started every 





25 April 0600 G'ff 
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Figure 1.9. jue -1.0 isopleth of programs for ESI, 25 April 0600 GMT. 
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50 mb and ranged from 800 to 250 mb. It was assumed, to maintain the 
same relative relationship of ESI among various programs, the SI 
obtained from Program 3 was divided by two, However, a weighted 
average may be more representative. 

The forecasting ability of Program 3 was tested using the same 
AVE IV data, The results are shown in Figures 1,7, 1.8, and 1.9. The 
line with the circled 3 is the *1,0 isopleth obtained <>om Program 3. 
The figures show that some improvement was made for all three time 
periods. Figure 1.7 shows that for 24 April 1600 GMT improvement was 
made In eastern and south-central United States, Stations 260, 240, 
34C, and 317 wpre excluded from the false alarm area, Figure 1,8 
shows that for 24 April 2400 GMT improvement was made in southeastern 
United States by eliminating Station 213, Figure 1,9 shows that for 
25 April 0600 GMT improvement was made in northern and southern United 
States, Stations 240, 226, 208, 425, 429, and 456 were excluded from 
the false alarm area, A specific list of stations gained or eliminated 
is ^iven in Table 1,1. The modification of the program involving the 
number and location of the upper level air layers had some beneficial 
effects , 

To determine the importance of layers located between the 
lifting condensation level and the 800-mb level to blocking potent-; ,1 
in double vortex thunderstorms, the upper level layers were shifted as 
a whole to begin at the lifting condensation level. The third modifi- 
cation (Program 4) used 12 layers and each layer was 150-mb-thick. The 
first layer begins at the lifting condensation level and a new layer 
starts every 50 mb. 
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Table 1,1, Comparison of Results Obtained from 
the Various Programs 


Stations Eliminated 


Stations Gained 


Program 2 
1800 GMT 24 
2400 GMT 24 
0600 GMT 25 

Program 3 
1800 GMT 24 
2400 GMT 24 
0600 GMT 25 

Program 4 
1800 GMT 24 
2400 GMT 24 
0600 GMT 25 


April 

April 

April 


Apri 1 

260 

April 

213 

April 

240 

429 

April 

April 

April 

226 


248 

340 

317 

226 

208 

425 

456 




456 


213 


340 425 
208 317 
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The forecasting ability of Program 4 was tested using the same 
AVE IV data. The results are shown in Figures 1.7, 1,8, and 1,9 by the 
line with the circled 4. Figure 1.7 shows that for 24 April 1800 GMT 
Stations 340 and 425 were added to the false alarm areas that were 
located in northeast and central United States. Figure 1,8 shows that 
for 24 April 2400 GMT Stations 226 and 456 were eliminated in northwest 
and southeast United States, Figure 1.9 shows that Stations 208 and 317 
were gained in the false alarm area in eastern United States. A 
specific list of stations gained and eliminated is given in Table 1.1. 

No significant increase or reduction in false alarm area was achieved 
by Program 4, The layers between the 800-mb level and lifting conden- 
sation level appear to have little effect on the forecasting ability of 
the ESI , 

Summary of Predictor Approach 

A program of Eagleman's model of a double vortex thunderstorm 
was modified three times. Each modification was to reduce the false 
alarm area and test the effect of forecasting ability of different 
layers in the double vortex thunderstorm. The forecasting ability of 
each modification was tested using the same AVE IV data. The results 
of each program were compared to the results of the previous program. 

The modification involving the low level layer showed little effect on 
forecasting ability. The modification involving additional upper level 
layers located at higher altitudes showed some beneficial effect on 
forecasting ability. The modification involving layers at lower 
altitudes (i.e., layers between the lifting condensation level and the 
800-mb level) showed little effect on forecasting ability. The results 


obtained during the sequence of tests that, modifications to Eagleman's 
tornado prediction program involving upper level layers at higher 
altitudes, had the greatest beneficial effect on forecasting ability. 
However, the results obtained are still deficient for operational use. 

As a result of the deficiencies of these various programs in 
predicting tornadic storms, the present investigator believes that an 
additional parameter should be incorporated into the index to account 
for other dynamical effects characteristic of double vortex thunder- 
storms. An identification of the parameter must come about through a 
better understanding of the fluid mechanics in double vortex thunder- 
storms. 


CHAPTER II 


FLUID DYNAMIC REPRESENTATION OF 
DOUBLE VORTEX THUNDERSTORMS 

The numerical tornado prediction scheme, explained and tested 
in the previous chapter, is based on the double vortex tnunderstorm. 

The present Investigator believes the scheme is incomplete and suggests 
that an additional dimensionless parameter be incorporated into the 
Energy-Shear Index (ESI), This parameter should represent the fluid 
mechanics of blocking in double vortex thunderstorms and be a measure 
of the shear instability associated with tornadic development. The 
parameter could be used to indicate a triggering mechanism in double 
vortex thunderstorms, A triggering mechanism initiates tornadic 
development and is needed to realize the blocking potential. The large 
"false alarm" areas could be eliminated if stations in the area fail to 
indicate a proper triggering mechanism. The purpose of this chapter is 
to demonstrate the feasibility of obtaining information on this 
parameter by wind tunnel investigation. 

Eagleman suggests that a thermal updraft of a double vortex 
thunderstorm blocks the upper level air layers as a cylindrical body [2], 
The flow around a cylindrical body separates on each side of this body 
creating a shear region behind it. The resulting wake of a thermal 
updraft is known as the dynamic updraft. It is in this dynamic updraft 
that tornadoes are believed to be created [2]. 

Turbulence levels in and around cumulonimbus clouds are 
extremely large [8], The blocking effect of a thermal updraft can be 
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simulated by a cylinder in the crossflow that has sufficiently large 
turbulence levels, The wake of the cylinder would be similar to 
dynamic updrafts in double vortex thunderstorms. The wake could be 
analyzed under laboratory controls and information on the above 
suggested parameter could be obtained. Turbulence levels required for 
complete similarity between atmospheric thermal updrafts in double 
vortex thunderstorms and laboratory flows about cylinders would be 
extremely difficult to create in a wind tunnel and are beyond the scope 
of this study. Instead, a wind tunnel investigation of a cylinder in 
laminar crossflow was conducted. The conditions created in the wind 
tunnel are represented by a Reynolds number (based on cylinder diameter, 
crossflow velocity, and kinematic viscosity) of 5,000. This is of 
equal value to a Reynolds number obtained from multiplying average 
storm diameter to average wind speed in storms and dividing by the 
amount of turbulent transport of vorticity measured by aircraft in and 
around cumulonimbus clouds. Alehsandrov [8] had made these aircraft 
measurements and obtained the value of: 

Re t = 2.0 x 10 3 to 1.0 x 10 4 , (2.1) 

2 

where a turbulent eddy viscosity of 20 to 100 rn /sec is used. 

Although the above Reynolds numbers are similar in value, it 
should be remembered that one is based on kinematic viscosity, while 
the other is based on eddy viscosity. Consequently, the present 
investigator does not suggest that any dynamical similarity exists 
between the atmosphere and the tunnel conditions. However, it was 
strongly felt that an investigation of the dynamic behavior of vortex 
shedding mechanics and/or shear wave instability behind cylinders in 
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the crossflow could shed light on the mechanism of tornadic generation 
behind thermal updrafts, The concept to be described, however, is the 
important issue. 


I. PREVIOUS INVESTIGATIONS 

Previous investigations of the phenomenon associated with the 
blocking effects of cylinders in a laminar crossflow have revealed an 
alternating periodic vortex street behind the cylinder. The first 
pictures of the alternating periodic vortex street behind a cylinder 
were published in 1902 [9], In 1911, von Karman presented his famous 
vortex street theory that explains a stable alternating row of vortices 
could exist in a street behind a cylinder [9]. All Reynolds numbers 
used in discussing the phenomenon observed and measured behind 
cylinders in wind or water tunnels are based on free-stream velocity, 
cylinder diameter, and kinematic viscosity. 

The Reynolds number range of periodic vortex shedding behind a 
cylinder is divided by Rosko [9] into two distinct subranges. The two 
ranges are the stable range and the irregular range. The stable range 
is composed of stable laminar vortices in a street extending for long 
distances behind the cylinder. The irregular range is composed of 
vortices in a street behind the cylinder where turbulent velocity 
fluctuations accompany the periodic formation of vortices. The 
diffusion of vorticity in the irregular range is by turbulent transport. 
This transport causes the wake to approach isotropic turbulence 40 to 
50 diameters downstream [9], 

At conditions represented by a Reynolds number range of 40 to 
150, the alternating periodic vortex street is in the stable laminar 
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range, and can be observed for long distances downstream, The laminar 
vortices are dissipated by molecular diffusion of vorticity only. The 
mechanism of laminar vortex production operates in a way that a 
Strouhal number increases with an increasing Reynolds number [9]. The 
Strouhal number represents the shedding frequency in non-dimensional 
form. Figure 2.1 depicts this relationship in a graph of Strouhal 
numbers versus Reynolds number [10], 

The Reynolds number range of 150 to 300 is a transition region 
bridging the stable range and the irregular range. The Reynolds number 

5 

range of greater than 300 to less than 3 x 10 represents the 
conditions of the irregular range. The vortex street for this range is 
associated with turbulent velocity fluctuations. The mechanism of 
turbulent vortex production operates shortly downstream of the cylinder 
and in a way that the Strouhal number remains clearly constant with 
increasing Reynolds number, as shown in Figure 2.1. 

In a portion of the irregular range of the von Karman vortex 
street a laminar shear region exists behind the cylinder. The end of 
the laminar shear region is unstable and alternately rolls up Wjlpn 
forming turbulent vortices that make up the von Karman vortex 
street [11], This laminar region begins at the separation points on 
the cylinder and ends at total breakup associated with transition. The 
total breakup, at conditions represented by a Reynolds number in the 
range of 3,500 to 8,500, occurs at (0.7) to (1.4) diameters downstream 
from the separation points on the cylinder [12]. 

Bloor [13] noted transition waves occurring in the laminar 
shear region behind the cylinder. The shear region consists of laminar 
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Figure 2-1. Strouhal number vs Reynolds number for 
a circular cylinder [9]. 
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separated boundary layers, These waves precede the total breakup to 
turbulence and are believed by Bloor [13] to be two-dimensional, 
Tollmien-Schlichting waves caused by the instability in the separated 
boundary layers. Bloor has measured the frequency of wave production 
at conditions represented by a Reynolds number in the range of 500 to 
50,000. 

Gerrard [14] photographed these waves at conditions represented 
by a Reynolds number up to 2,000 by a flow visualization technique. 
Figure 2.2 is a drawing of his results, The linear portion of the 
transition waves can be seen to exist between three and somewhere near 
five radii downstream of the cylinder's center, Beyond this, the 
breakup to turbulence and the periodic formation of vortices that make 
up the von Karman vortex street exist. 

Gerrard's flow visualization technique allowed transition waves 
to be photographed at a Reynolds number up to 2,000, To investigate 
transition waves by a flow visualization technique at conditions 
represented by a Reynolds number greater than 2,000, a wind tunnel was 
constructed. 


II. CONSTRUCTION METHOD 

The wind tunnel was constructed in a 6-ft by 30- ft by 8-ft 
trailer. One-fourth-inch-thick corrugated center cardboard formed the 
wind tunnel test section and the inlet nozzle, Three-inch-wide masking 
tape was used to fasten the cardboard edges together. Two layers of 
1/16-in. mesh screen covered the opening of the inlet nozzle. The wind 
tunnel has a contraction ratio of 5.5, and its dimensions are shown in 
Figure 2,3, 
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A 1/3-hp electric motor was used to turn a squirrel -cage fan at 
the opposite end from the Inlet nozzle, The fan sucked air through the 
tunnel and expelled it to the outside, Materials were needed to insure 
the straightness of airflow. Incorporated into the tunnel design were 
eight 3- in, by 6- in, aluminum cans, The ends of the cans were removed 
with a hacksaw so that they resembled a tubular structure. The cans 
were positioned immediately upstream of the squirrel-cage fan. In 
order to produce a honeycomb-like structure, experience Indicated that 
straws produce the closest approximation. Therefore, upstream of the 
aluminum cans was composed of 1/4-in. -diam plastic soda straws, 1.5 in. 
in length. Held into place by two 1/16-in. mesh screens on each side, 
the straws were stacked to resemble a honeycomb. The test section is 
painted black with white lines in a square grid. 

In order to view the test section, clear glass windows were 
incorporated in the tunnel design. Windows (21 in. by 16 in.) were 
positioned on the side and top of the test section. It was desirable 
to observe the test section from both the side and top views simul- 
taneously. Consequently, a 24-in. by 18-in. mirror was placed over the 
glass on the top of the test section at a 45-deg angle from the 
horizontal plane. It provided a simultaneous view from the side and 
an inverted top view. 

A titanium tetrachloride smoke generation system was developed. 
The system combines the use of a smoke gun and a fish tank air pump. 
White, neutrally buoyant smoke was pumped from the smoke gun through a 
3/16-in. -diam plastic tube. At the end of the plastic tube is the 
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Injection nozzle. The nozzle was made from a 0,02-in. -ID stainless 
steel tube. 

Two slide projectors were used to Illuminate the test section. 
An Olympus OM-1 35-mm single-lens reflex camera was positioned at the 
side of the test section such that the field of view included the 
simultaneous inverted top and side view. 

Flow visualization techniques require wind tunnels with 
extremely low turbulence levels. Extreme care was taken to reduce 
these turbulence levels. The tunnel is located in one end of a 6-ft by 
30-ft by 8-ft trailer, A wall divides the trailer into two sections. 
The section containing the tunnel consists of approximately one-third 
of the trailer's volume and is sealed off from the outside air. The 
dividing wall contains a door with a 3-ft by 3-ft section removed. Two 
1/16-in. mesh screens, 1 in. apart, cover the removed section. This 
allows low turbulence air to enter the portion of the trailer that 
contains the tunnel, A door that opens to the outside air is located 
in the remaining two-thirds of the trailer. Since the air that flows 
through the tunnel is expelled to the outside during testing, this door 
was opened slightly to allow air to enter the trailer. 

Initially, a velocity of 2.8 ft/sec (i0.2 ft/sec) was measured 
by a constant temperature hot-wire anemometer throughout the test 
section. The smoke was injected in the flow at the center of the 
beginning of contraction of the tunnel's inlet nozzle. Figure 2.4 is a 
photograph of smoke along the central streamline of the test section. 
The shutter speed is 1/60 sec. The straightness and sharp definition 
of the smoke line is evidence that the test section is laminar and in a 
steady state. 
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a. Top View 



b, Side View 


Figure 2.4. Photograph of center streamline of tunnel 
(shutter speed 1/60 sec). 
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A test stand that allows the testing of an air jet and a 
cylinder in crossflow was built. This stand and its dimensions are 
shown in Figure 2,5. The horizontal sheet of the test stand is 
1/32-in. galvanized sheet metal. The sheet completely spanned the test 
section to divide it horizontally into an upper and a lower half. A 
12-in. square section of the sheet metal directly behind the cylinder 
tube was removed and replaced with 1 /8-in. -thick clear plexiglass in 
order that the lower half of the test section could be viewed from 
above. The sheet metal was held in the center of the test section by 
clear plexiglass sides. 

The test stand is designed such that a cylinder is positioned 
in the lower half of the test section. The cylinder is a 3.5-1n.-diatn 
hollow cardboard tube. The tube has a 1/4-in. wall thickness and is 
assumed to be round to within 0.02 in. The surface of the tube is as 
smooth as paper and is painted black. Plastic straws (1 .5-in. -long by 
V4-in.-diam) are placed in the bottom of the tube to resemble honey- 
comb. When the test stand is mounted in the test section, an air jet 
of varying speeds may be expelled through the tube into the upper half 
of the test section. When the tunnel test stand is installed in the 
test section, the flow in the test section is measured by a constant 
temperature hot-wire anemometer throughout the section to be 
3.0 ft/sec ±0.2 ft/sec. 

The Reynolds number representing the conditions of flow about 
the cylinder at the test conditions is computed by: 


Re 


VxD 


5,468 , 


v 


( 2 . 2 ) 


Top View 



^isure 2.5. Sketch of test stand. 
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where 

V - 3.0 ft/sec, D * 3,5 In., and v « 160 x 10"^ ft^/sec . 

The value of the kinematic viscosity is obtained from Schlichting [15], 
The bottom end of the cylinder rests on the bottom of the test section. 
The boundary layer is assumed to remain small in the favorable pressure 
gradient created by the contraction of the tunnel's inlet nozzle. The 
cylinder is located 5 in. ■‘wnstream of the end of the contraction of 
the tunnel inlet nozzle. 

The horizontal section of the tunnel test stand holds the top 
of the cylinder in place and extends 5 in. upstream of the cylinder. 

The height of the laminar boundary layer on the test stand and tunnel 
wall at the location of the cylinder is calculated from [16]; 

6 = t = 0- 3 in- (2.3) 

V ] '< c 
00 

vX 

The grid painted on the wall of the test section is one radius 
of the cylindrical tube that was incorporated in the design of the test 
stand. The distance X is measured in the downstream direction from the 
center of the cylinder, while V is measured in the cross stream 
direction from the center of the cylinder. The origin of the grid is 
the center of the cylinder. 


III. TEST PROCEDURE 


To visually display the waves in the shear region behind the 
cylinder, smoke was introduced in the flow at the beginning of 
contraction of the tunnel's inlet nozzle. This was done so that the 
smoke line encountered the cylinder just off center of the front 
stagnation line on the cylinder. A photograph taken with a very short 
exposure time will display a smoke-filled laminar streakline as if it 
were stopped in time. Figure 2.6a is a photograph taken with the 
camera's shutter speed set at 1/60 sec. The field of view of the 
photograph extends over the range of X = 1 to 7.5 radii and of 
V = 4.0 to -2.0 radii . 

Due to its turbulent nature, evidence of the von Karman vortex 
production rate is not shown on film. However, the turbulent region 
past X = 3 . 0 radii was visually observed to oscillate near 22 times 
per 10 sec. This corresponds to a Strouhal number of 0,21, which is 
comparable to previous results by Roshko [9]. 

The prints shown in Figures 2.6a, b, c, and d were reproduced 
from slides and have in the process lost some of their sharpness. 
However, the definite sharpness of the smoke line indicates the smoke 
streakline up to 2.5 radii downstream of the cylinder's center is 
laminar. There is no problem in distinguishing the difference in 
sharpness levels before 2.5 radii and after 2.5 radii downstream of the 
cylinder's center. 

The frequency of wave production was estimated by taking photo- 
graphs at a variety of shutter speeds. Figures 2.6b, c, and d are 
photographs of the same region taken at camera shutter speeds of 1/30, 
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I. Shutter Speed 1/8 sec 
Figure 2.6. (Continued) 
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1/15, and 1/3 sec, respectively. Figure 2,6c shows movement of the 
wave on the smoke streakline making it less distinguishable than 
Figure 2.6b, while the wave in Figure 2.6d is completely indistin- 
guishable, This indicates a wave production rate of greater than 
8/sec and less than 15/sec. The ratio of wave production to von Karman 
vortex production is between 3.6 and 6.8. This estimated ratio 
compares well with measurements by Bloor [13]. The test results ratio 
and Bloor's measurements and ratio are shown in Figure 2,7. 

The photographed waves in the laminar smoke streakline before 
3.0 radii downstream of the cylinder's center, at conditions repre- 
sented by a Reynolds number near 5,000, appear to be remarkably similar 
to the linear waves before 3.0 radii photographed by Gerrard. At 
conditions represented by a Reynolds number of 5,000, the breakup of 
turbulence appears to occur near 3,0 radii downstream of the cylinder's 
center. In Gerrard' s photographs, at conditions represented by a 
Reynolds number of 2,000, the breakup to turbulence does not appear to 
occur until near 5.0 radii downstream of the cylinder's center. In the 
additional 2.0 radii that the laminar shear layers exist, at conditions 

represented by a Reynolds number of 2,000, the transition waves become 
non-linear and can be seen to have rolled up. At conditions repre- 
sented by a Reynolds number near 5,000, the non-linear region, where 
transition waves roll up, could not be seen by the described flow 
visualization technique. 

To see if the same trends occur in the wake behind an air jet 
directed normal to the flow, the same procedure was attempted. The 
wake behind the air jet was turbulent, however, and no smoke lines 
could be photographed. 
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Figure 2,7, 


Ratio of transition wake production to von Karman 
vortex production vs Reynolds number [13]. 
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IV. CONCLUSION 

Although no relation can be drawn about the wake behind the 
tested cylinder and double vortex thunderstorms, the feasibility of 
obtaining information from a wind tunnel investigation was demonstrated. 
Further research is needed to develop a means of creating levels of 
turbulence in wind tunnels sufficiently large to guarantee similarity 
with natural conditions. If such a means can be created, then, the 
wakes of cylinders may be studied under laboratory controlled 
conditions. Information about the fluid mechanics in double vortex 
thunderstorms and instability associated with tornadic development can 
be obtained and incorporated into a dimensionless parameter. Such a 
dimensionless parameter could be incorporated into the Energy-Shear 
Index and improve the Index's ability to predict tornadic development. 
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PROGRAM 4 




«uqd t i he^. pages »?9 ?hol iut 

C THIS PROGRAM CALCULATES THE CNCRGY-SHEAR INDEX FKDN RAUIHGOHDG DATA 
C DATA IB READ FRUM CARDS WITH THE FULLOUING FORhA. 

C COES' )-7 GEOPFNTIAL HP l OUT i ODES 0-14 PRESSURE III MP> 15-00 TEMPERATURE C 
C COLB 21-26 DEWPOINT 1EMI'. C. 2V-32 WIND H ! IVfCT ION . 3 3 - 3 V WHIP SPEED M/BEC 
C 30-43 EASTWARD COMPONENT OP UIHPi 44»49 HURT HOARD t OMPOHENT OK WIND 
B 50-55 POTENTIAL TEMPERA TURK i S6-A1 EOUIVAUIIT POTENTIAL TEMPERATURE 
C COL 62-64 MIXIHO RATIO 
C COL 67-BO ARE USED POR IDENTIFICATION 
C 

C EACH DATA SET IB PRECEDED PT AH ID CARD WITH THE FORMAT 
C COL t-5 STATION NUMBER. 6-33 STATION NAME. 34-35 DAY . 36-45 MONTH. 

C 46-4? YEAR. 40-52 TIhE OF DAY 
C 
C 
C 

C DTEHP» DEWPOINT TEMPERATURE I TEMP- TEMPERATURE (DEO C>( ET “ TOTAL ENEROT OF 
C AH AIR PARCEL I WB " WIND SPEED AS MEASURED FROM THE SOUNDING) WD « WIND 
C DIRECTION A8 MEASURED) IOTA « STATION NUMBER I 0PH1 b OEOPOTCNT1AE HCIONT IN 

C METERS) AMIXRA* MIXING RATIO <OM/KOM)l El » ENERGY IMD) X 

DIMENSION DTEMP( 3) » DATG(&) i PRESS (3) iOPIIT(3)» TEMPO) »KT<2) »UB(41 1 
DIMENSION WIt(4i ) » DSRHiD<12.200) , BPBHI D < 1 2 . 200 > i RIJS < 40) , RWD ( 401 
DIMENSION CUS <40 1 .MfSB25< 12. 200 1 .MSD25! 12.200) .68253( 12.200) 

DIMENSION 6)12502.200) .ISAVEI 12). DCU25U2, 200) 

DIMENSION SRW025 0 2;200)iSRWD2U0 2.200) 

DIMENSION HEIGHT (25) »PRESH<25) 

INTEGER HOUR 
C 

DATA AUEFT/4HLEFT? 

DATA RIGHT/4HR0HT/ 

C 

COMMON T CUP * DTEMP. WD. WS. HOUR. NDAYo PRESS. CPHT 
COMMON IOTA 
READ IN A DATA BET 

100 READ (50102) IOTA. DATE 
IF OSTA .EO, 9990?) STOP 

• READ <5.004) QPHT(l). PRESS U). TEMP <1>. DTEMP ( 1 ) t WD<1 ), WS< 

1 1 1 .HOUR. HDAY 
DO 10S LDP1“ 2.7 

READ (5. 1103 MIEIOHT (EOf'l >» PRESH (LOP) ) . UP ( L0P1 > » US(LOPl) 

105 CONTINUE 

READ (50104) (IPHT<2). PRESS <2). TEMP <21. DTEMP 1 2 ) . WD<0 ). WS< 

1 0) 

HEI0HT<0)nG PHT < 2 ) 

PRESH<Q)»PRES8<2> 

DO 110 LOPS'. 9.21 

READ <5. 0 03) HEIGHT ( LOF'2 ) t PRESH (L0P2) .WD(L0P2>i WS(L0P2) 

OO CCHIT I HUE 
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READ <5,1104) 01 H T (3 ) i PUEM (3), TEMP (3), OTCHfi J> » Ut«< V? > » VP< 

J 22 ) ’ 

PCM 15 LOP3* 23 Mi 
READ < 0 » 1 1 0 A > HU < LOr-3 > r H5H.0P3) 
no CONTI HOC 
1102 FORMAT < 1C . 2UX ,t5A4 > 

1103 FORMAT <2I’7 . 1 , 1 C'X > I 4 (< , 1 » f 0 , 1 > 

1104 FORMAT <?F? . 1 , 3F6 , 1 i F5, 1 , 35X » 10 , 13 ) 

1106 FORMAT (26X, P6 , 1 » Ft,, 1 ) 

na*tTi:MP(i)-itiTt;HP<i))/o.tiooo,ioPHr<i> 

C IF PATA IB HISSING, PRINT MESSAGE 

IF < 0Pliri2) ,EU, 99.9 .OR, QPItTf 3) ,E0, 99,9) 00 TO 299 

IT ( TCMP(2) ,EM, 99,9 , OR , TEMPO) ,F0, 99,9) 00 TO 299 

IF { DT£HP<2 > ,E0, 99.9 , OR , UTEMP! 3> ,E0. 99,9) 00 TO 399 

C 

c 

C CALCULATE THE ENERGY INDEX 

C 

BO 240 If, 2, 3 

TCMf'A * PTEMP 11) +273,16 

TEMPO * TEMP (I) +273.16 

IF (TEHPA . LT < 273,16) GO TO 220 

EXP = <17,2693002 * < TEMPA-273 . 1 6 ) >/< TEHPA-35, 06) 

GO TO 230 

220 EXP * ( 21 , 074000 + < TEMPA-273 , 1 6) )/ ( TEMP A “7 , 66 ) 

230 E « 6, 10704(3, 710284+EXP) 

AMIXR’A « <622 , +E )/< PRESS < I )-E ) 

J*I-1 

ET ( J) * , 24Y(TEMPB +2 , 54AK1XRA+CP1IT < I ) /.1 00, > 

240 CONTINUE 
. El* ET< 2 ) "ET < 1 ) 

00 TO 300 

299 UNITE (4,1012) HATE, I8TA 

1012 FORMAT ( 42UTME ENERGY INDEX CANNOT PE CALCULATED FOR ,3A4,1X»I5> 

00 TO 100 
C 

C**+l!+.t.!M*W****4**lM ***m,»**»4*********»»»*4**Y* + *I’*4***M'***l.*****# 

C 

C CALCULATE THE MEAH HIND VELOCITY OF THE • ’ OUD LAYER 
C USING THE UIHDB AT THF 830,700,000,1 300- MB LEVELS 
C 

300 WRITE (6,1122) DATE, IBTA 

1122 FORMAT (7////// , 1 7X, 5A4 , IX 1 15, 40X, 10IISHEAR INDEX VI H0D./7) 

C THIS LOOP TOTALS THE NUMBER OF 50-MB LEVELS THAI HAVE NO DATA AVAILABLE (XCI) 
C AND THE NUMBER OF LAYF.RB (USED TO CALCULATE THE AVERAGE) WHICH HAVE HO DATA 
C AVAILABLE! IFOUR) 
tC I * A 

r : uuR * o 
DO 304 1*1,20 
DH 

rr (I , NE. 1) IH»2*I-2 
IF (USt IM) ,NE. 99,9) GO TO 304 
JCI« lCl+1 

IF (IM.Cfl.O , OR , IM.E0.14 ,0ft, IM,E0,22 ,0R. IM.EQ.30) CO TO 302 
GO TO 304 

302 IFOUR * I FOUR 4 1 
304 COHTINUE 
306 M * 20-rCI 
XUINB =0. 

YU I HU * 0. 

C THIS LOOP HILL BUM THE HINDS AT THE 050,700,500,1 300-HB LEVELS 
ti DIVIDE DY rur. HUMBER OF LAYFRO TO 0DTA1N I HE AVERAGE , THEN FIND THE 
C MAGNITUDE (THUS > AIIU DIPECT10H ( TMWD) IK DEGREES 
UP 312 J*'J ,20 
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IF (USMM) .EClt 99,9 > 00 TO 312 

IF (JM.FO.D i OR i JM.Efl.14 .OR, JH.EU.22 ,015. JM.C0.30) CO TO 310 
00 TO 312 

310 RADWS p UIK JH>/5?,296 

XU 1110 XU1ND f UO(JM)lCOO(R.niUO> 

YW1ND p YU I HO 4 UB( JMJIGIMIRADWS) 

312 CONTINUF 

320 AVXU n XWlND/<4 .-FLOAT < lF(HIR) ) 

AVYU * YUIND/( 4 . -FLOAT! H OUR ) ) 

THUS® SflRT(AVXU+4? FAVYWm') 

00 1= AVYU/AVXW 

TMUD « ATAM100MS7 ,296 

C FIND THC REAL VALUE or TMUO OY SELECTIMO THE PROPER OUADRANT 
IF< AVXU , OF . 0 . ) 00 TO 350 
THUD r. THWD +100, 

350 IF (THUD ,LT. 0.) THUD - THWHF360. 

IF ( THUD ,0T, 360.) THUO '• THUD-340. 

C INITIALIZE VALUE OF ISAVE AHD IHDSHR 
00 365 M»i,ll 
ISAVE(H) *1 
365 CONTINUE. 

IHDOIIR «• 0 
no 699 m« 1,2 

DO 699 K « 1 ,61 ,5 
DO 499 I. » 1,31,5 

C THE VALUE OF III DETERMINES UHETDFR DEVIATION JO TO THE RIGHT OR LEFT OF THE 
C AVERAGE CLOUD LAYER UIND/K IS THE ANGLE OF DEVIATION. L+SO IS THE 
C PERCENTAGE OF AVERAGE WIND SPEED 
C 
C 

c 

C CALCULATE RCLATIVE WIND VEIOCITIES DASED ON DEVIATIONS 
C OF STORM DIRECTIONS FROM THE HCAN WIND .VELOCITY 
C 

IHUNPRn S04L-1 
IDEVIT * K-l 

C 05* TEST STORM SPEED! SO* TEST STORM DIRECTION 
SS » TMUS IFLOAT(50FL-l)/100. 

IP (III.EO.n DIR” RIOHT 

IF (III.EO, 2) D IR“ ALEFT 

IF ( niR.Efl, ALEFT) IDEVIT= IICVIT+5 

IF (III, £0.2) 00 TO 401 

SD=> THUD 4 FLOAT (l> >**1 . 

00 TO 402 

401 Bit* THUD-FLOAT (K)-4* 

402 DO 430 LL= 1 , 20 

LLM = LL 

IF (LL.NC.l)LLH*2*LL-2 

C CAMMA ” DIFFERENCE DETWEEH TEST STORM DIRECTION AND DIRECTION OF MEASURED 
C UINDI RADOA «* OAHMA IN RADIANS 
OAHMA “SD-Wn ( LLM ) 

RADOA" GAMMA/57,296 

C RWS « RELATIVE WIND SPEED DETWEEN THE STORM AND THE MEASURED WIND 
RUSULH) = S0RT(WS(i.LM)**24SS*t2-2.0*US(LLH)*SS|C0S<RADGA) ) 

SINGS " ( 05-fSIN ( RADGA > ) /RUS (LLM ) 

C. THIS SEGMCNT PREVENTS THE CENTRA! ION OF SINES OF VALUE LARDER THAN 1.0 IN 
0 ADSOI.UTE VALUE DUE TO ROUND OFF ERRORS IN THE COMPUTER 
NS IN -"0 

IT (SINES, LT. 0,0) NSIII«1 

IP (ADS < GT NHS> • GE , , VV9V99R .AND. ADS ( BINOD , LT. 1 . 0001 > SINSS=1,0 
IF (SlNSS.bT, 1. ! GO TO 1224 

IF (NS1N . E(t > 1 , AND i DINS', ,E(I, 1,0) SINGS « -1,0 
IP ( ( 1 , O-ADS ( D I N S3 ) T*2 , > ,0T, 0.0) 00 TO 405 
WRITE! 6, 5302) II I , K, L, LLM , S1NSS, COSOS 
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C0SSP»0,0 
DO 10 410 

405 C0C00 “SORT < 1 • 0— ADS <SINSS>**2 1 0) 

5302 FORhAT <4I4,2E14,4) 

410 IF (COCOS ,LE. ,0000011 00 TO 415 
TANSO” SINGD/C0508 
SETA « ATAN< TANSS 1 *57 , 294 
DO TO 420 

415 WRITE (6, 53021 II I > K > L , LLMr SINSB, COSSS 
IF (SINSG - 0,0) 4U,4U,41? 

414 SETA « 270, 

00 TO 420 
41? SETA " 90, 

420 RWtULLHl “WO (LLM1-BETA 

IF (RUD < LLH) . LT . 0, 1 RUD(LLM)= RUIKLLMH340, 

IF (RUD(LLH>.GT. 340,1 RUP(LLH>“ RU1XLLM1-340, 

' 430 continue 

C N IS THE NUMBER OF 50-HD LEVELS THAT HAVE DATA 

C CALCULATE THE AVERAGE RELATIVE WIND IN THE SURFACE-850HB LAYERI 
DO 435 1CI«2, 10,2 
100*10-1 
IIGuIQOf? 

1110*100+29 

IF (HEIGHT < 10) . OT, HOT GO TO 437 
435 CONTINUE 
437 ILCL«IOO?2+2 

ini * o 

D0440 H*1 t I LCL 
MM“H 

IF (H.HE. 1) HH*2*H-2 
IF (WS(HH).NE, 99,9) 00 TO 440 
IPI * IDI+1 
440 CONTINUE 
XSUH * 0, 

YSUH “ 0, 

NN=ILCL-IDI 
DO 450 1 1* 1 , ILCL 
I IH=I I 

IF ( II.NE.l) I IN* 2*11-2 

IF ( UBUIMl.EQ. 99,9) GO TO 450 

RRUD * RUD ( IIM) 757.294 

XSUH= XSUH+RUS< IIM1*C0S (RRUD) 

YSUH *'YSUH+RUS(IIH)*SIN<RRUD> 

450 CONTINUE 
C 

C********* ************************************** *********************** 

c 

C CALCULATE LAYER AVERAGES 
♦ C 

C AVRUS AND AVRUD ARE THE AVERAGE RELATIVE WIND SPEED AND DIRECTION 
C RESPECTIVELY OF THIS I.AYER 

C THIS PROGRAM SEGMENT CALCULATES THE SURFACE- 850-KB AVERAGE KIND SPEED AND 
C DIRECTION 

• AVXSUH “ XSUM/FLOAT <NN> 

AVYSUM = YSUM/FLOAT ( HH 1 

AVRUS “ SORT ( AVXSUH**2+AVYSUM**2) 

D “ AVYSUM/AVXSUM 
AVRUD r AT AN (D 1*57, 294 
IF (AVXSUH. OE. 0,1 GO TO 501 
AVRUD * AVRUD+100 , 

501 IF ( AVRUD ,LT , 0,1 AVRUD * AVRUD+340. 

IF (AVRUD ,GT. 340,1 AVRUD “ AVRUD-340.0 
CUD = AVRUD +180, 

IF ( CUD.CT. 360,1 CUD <= CUD-340. 

DO 570 KK“I IQ r 1110,2 
KKJ '0 
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XRU "0, 

YRU * 0. 

KKI *= RK-6 
NKtfnKK 

DO 520 KKL H KKI.KKKi? 

IP (US(KKL) >F.0*99,9) BO TO $10 

RRRWP «RWP < KKL ) /67 >296 

XRU « XRU+RUS<KKLMC03(1«RRUD) 

VRU u VRM + RUS < KKL > * B I N t RKRUD ) 

GO TO $00 
510 KKJ= KKJ+i 
620 CONTI HUE 
525 IF (KKJ.EC1.4) 00 TO 550 

XAVRU * XRU/ ( “1 > -FLOAT < KKJ > 1 
YrtVRM ■ YRU/(4,-FL0AT<KKJ>) 

' XAVRUS * SORT (XAVRUM'Sf YAVRU**2) 

F.» YAVRU/XAVRU 
XAVRUD a ATAN(E>*57,296 
IF ( XAVRU, 0E,0.) O0 TO 630 
XAVRUD » XAVRUD+ 100 i 

530 IF ( XAVRUD . t,T 1 0 . ) XAVRUD * XAVRUP+3A0, 

IF < XAVRUD • OT , 3A0.) XAVRUD " '(AVRUD-3A0, 

CDE61 “ <CWD-XAVRUP> 

IF (CDEOl , OT , 90 . > 00 TO 560 
CHEG2 * CPE61 + 90. 

IF (CPEG2 i l.T .0 1 ) 00 TO SAO 
540 CHEG = CDEOl/57,296 

CU8 ( KK ) «XAVRUStCCIS (CDEO) 

00 TO 570 
630 CUS (KK ) =-1, 

CO TO 570 

560 XAVRUD a XAVRUD + 360, 

CBE01 « XAVRWD-CUD 
XAVRUD « XAVRUD-360. 

IF (CPE01.LT, 90.) 00 TO 540 
CUP r. CUD +360. 

CDEOl a CUD -XAVRUD 
CUD « CUP-3A0, 

IF (CDF.Ol ,LT. 90,) 00 TO 540 
CUS <KK>» -1. 

570 CONTINUE 
YhIP =0. 

XHID so, 

JJJJoO 

C THIS F'ROGRAH SEGHENT CALCULATES THE 300- TO 500-HD AVERAGE WIND SPEED AND 
C DIRECTION 

DO 560 JJJ«20r30,2 
IF (UBUJ.I) iEp. 99.9) GO TO 575 
RRRVID = RUD( JJJ)/57.29A 
XHID = XHID +RUS< JJJ)*COS (RRRIID) 

YHID s Ytim +RUS( JJJ)*S1N(RRRUD> 

00 TO SCO 

575 JJJJsJJJJ+1 
500 CONTINUE 

IF (JJJ.I.EO, 6) 00 TO S90 
XAVHID = XHIM/lA.-FLOAT(JJJJM 
YAVM1D = THlD/(6.-FLGAT(JJJJ>) 

• AVMOPD = SORT (XAVHID **2 + YAVHIB**2> 

0 s YAVHID/XAVHID 
AVHDIR = ATAtl ( 0 ) *57 . 29A 
IF ( XAVHID IliC . 0 . ) 00 TO 583 
AVHDIR s AVHDIR +100, 

505 IF ( AVHDIR . L T i 0 , ) AVHDIR = AVHDIR +3A0. 

IF < AVHDIR .OT, 360.) AVHDIR » AVHDIR-360. 

00 TO 600 
590 ftVHGPD a -1, 
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AVHU1R » -1, 

C 

c 

C CALCULATE SHEAR INDEX 
C 

600 JC25C » 0 
r r. 1. 

IXC2SC b 0 

DO 665 JJK«mi*III0i2 

ir <cus < j,ik> • l r 1 0 i > go to 6p d 

DIFF b POfi(.IJK) -AVRIIS 
A0C2 » .25TAVNU0 

IF <ABS<PIFF> .UEi ABC2> 00 TO 635 
625 F b 0* 

00 TO 665 

635 IF < F ■ EO i 0.) 00 TO 660 
TC25C * IC20C +1 
00 TO 650 
660 IC2SC *“ I 

650 IF < IC25C iGT . IXC25C) IXC25C » IC2SC 
F ■ l . 

00 TO 665 
665 CONTINUE 

IF < IXC25C » GT i INPSIIR) INOSHR - IXC2SC 
IF < IXC25C .EU. 0> GO TO 670 
00 TO 690 

670 IF (III. EG, I .AND. K.EO.l .AND. L.EO.i) 00 TO 690 

00 TO 6?9 
600 I ZERO « 1 

HH * ISAUEnXC25C> 

I1BS2S ( IXC25C »HH > « IHONPR 
MSD25 < IXC25C.I1H) « IDEVIT 
DEV23 < IXC25C . HH ) » DIR 
SRUS25 OXCQDC.HM) « AVRWS 
BKU025 (IXC25C.HM) » AVRUD 
SS25< I XC25C. HH ) =■ OF. 

S025 < IXC25C . HH > “SD 
DIRMID (IXC25C.HH> * AVHDIR 
SPOMIDUXCnSCfllH) - AVMSPD 
I8AVEUXC250 » ISAVE( IXC25C) +1 
GO TO 699 
690 IZERO =0 

ISAVE8 b IHUNFR 

1 SAVED « IDEVIT 
SAVEDV «= DIR 
GAVEH8 = AVRUB 
SAVEUtT b AVRUD 
^AVER BPg. 

SAVER" til 
SAVEHD b AVHDIR 
GAVEMS b AVHSPD 

699 CONTINUE 
C 

m*mr****t»***mm*m****»W *|.*******t*-*M. 
C CALCULATE ENERGY-SHEAR INDEX 
C 
C 

£BI b 6. -FLOAT ( INPSHR) / 2 .42 .Tfc'l 
IF ( I7ER0 . EG . 0) CO TO 920 
IlO 07.0 I Bill! 

if < isAvem.co, t> go to 9to 
MHN* ISAVL t I > -1 
C 

C OUTPUT RESULTS 
C 
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WRITE ( 4, 1 100 XIEIOHT ( 1UU > rPRESH U UCI ) , l(E10HT< 0) 

1100 FORMAT < ' ' , 'HEIGHT OF C P' ,F 1 0 , 2, tiX, 'PRESSURE OF 11 D',F10.2,5X, 
l'OSO MR HEIGHT' ,F10. 2) 

WRITE (4,1110) I 

1110 TORHAT ( // 1 33X , 44HTHE FOI.EOUING STORM SPEEDS AND DIRECTIONS CAVE A 
1CHEAR INDEX OF , It,/,') 

WRITE (4,1112) 

1112 FORMAT (SOX, 'STORM SPEED' , 2X» ISIISTORM DIRECTION, 3X, 3IIDCV, 4X, 

1 lAHBURF-USO AVCRAOE,CX, 1SH550-300 AVERAGE ) 

WRITE <4,1114) <HSS23U, J) , SB2SU , J) ,h9D2S< I , J> ,8D2S< I , J> , DEV20 
Kl,J>,liRWS2liU,J),SRUD2S(I,J) (SPriHIDa, J) ,D1RM1D<1, J),J»1,MMH> 

1114 FORMAT <28X> I?, 1" ,F5, 2, AH KPS, 3X, 12, IX, F6t 2, AH DEO , 3X, AS, 3X, 

1F5 1 2 , AM NFS, IX, FA, 2, 4|l DEC , 3X , F4 < 2 , AH KPS , 1 X, FA, 2, All DEO) 

GO TO 820 

910 URITF ( A , I 1 1 A ) I 

11 1A FORMAT (// , 31X, ' THERE UERE NO STORK SPEEDS AND DIRECTIONS THAT 
10AVE A SMEAR INDEX OF', 11) 

020 CONTINUE 
GO TO 930 
920 WRITE <A»U18> 

1110 FORMAT ( // » 4BX»36HTIIIS G0UND1NQ HAD A SHEAR INDEX OF 0,//> 

* WRITE (A, 1112) 

WRITE (A, 111 4) ISAVE3,SAVES,ISAVED,SAVED,QAVEDV,BAVEWS,SAVCWD, 

1 SAVEMD , SAVEHS 

930 WRITE < A , 1120) El , IND5HR,EDI 

1120 FORMAT ( // , 47X, 5HE1 “ ,FA.2,7H, SI » ,I1,12H> AND ESI «= , FA, 2/7) 

00 TO 122S 

1224 WRITE (4,4220) 1 1 1 , K , E , LUH , THUD , THUS , SI NSS 

4220 FORMAT <4IA,3E1S.7> 

1225 GO TO 100 
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GENERAL WEATHER CONDITIONS 

A tornado watch was issued for an area extending from central 
Tennessee into West Virginia at 1800 GMT 24 April. Severe thunder- 
storms and tornados occurred in northwest Tennessee, as a squall line 
moved into that region from Missouri, At this time, radiosonde 
stations in Nashville, Tennessee, and Marshall Space Flight Center, 
Huntsville, Alabama, recorded the highest wind speeds in the network 
at the 850-mb level. 

By 2100 GMT, a strong squall line had formed over central 
Tennessee into western West Virginia. The cyclone in central Kansas 
continued to intensity and moved southeastward into northern Oklahoma 
with severe thunderstorms actively developing rapidly along the cold 
front. At the 200-mb level, high-speed winds extended from eastern 
Texas, across northern Louisiana, into Tennessee, Virginia, and West 
Virginia. Highest speeds were recorded at Nashville and Huntington. 

In the upper atmosphere the two jets had further diverged, moving into 
the northern tier of states and into central Texas and southern 
Louisiana, respectively. 

Activity over the Appalachians had weakened by 0000 GMT 
25 April. The severe squall and tornado watch area was extended into 
southwest Missouri as the frontal system moved slowly eastward. Severe 
thunderstorms and tornados occurred over Oklahoma and southwestern 
Missouri . 
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By 0600 GMT, the cyclone had moved eastward into Kentucky, as 
the severe squall line grew in both intensity and length. Severe 
thunderstorms and tornados were occurring in eastern Arkansas and 
western Tennessee. Warm air was flowing northeastward from the Gulf 
in this area at the 850-mb level. At the 200-mb level, flow was 
strongly turbulent over the west Tennessee-west Kentucky region. 


APPENDIX C 


UNUSUAL SEVERE WEATHER PHENOMENA 
DURING EXPERIMENTS [7] 

1. April 24 2400 GMT, Grove County , Kansas: hall 2-1n.-diam 
(tennis ball), $50,000 to $500,000, 

2. April 24 2315 GMT, Wewoka, Seminole County, Oklahoma: hail 
3-in,-diam (baseball), $50,000 to $500,000, thunderstorm 

t = northeastern, 3 mi/8 mi. 

3. April 24 2325 GMT, Stilwell, Adair County, Oklahoma: hail 20 min, 
$50,000 to $500,000, 

4. April 24 2400 GMT, Cotton County, Oklahoma: hail l-in.-diam, 
2-ft-deep, $50,000 to $500,000, thunderstorm 7 = east, 3 mi/20 mi. 

5. April 24 2400 GMT, Craig and Ottawa Counties, Oklahoma: 
tornado $ = east, 5 mi southwest Miami, $62,000 to $500,000, 

6. April 24 1100 GMT, Nixa, Christian County, Missouri: thunderstorm 
? = southeast, two tornados (small), one east HW 65 and north 

Rt. 66, second Twin Acres. 

7. April 25 0040 GMT, Newton County, Missouri: tornado V = east, 

Rt. 60, south Neosho, $105,000, 400 yd +-+ 100 yd/9 mi. 

8. April 24 2400 GMT, Pettis County, Missouri: hail l-in.-diam, 

2- to 3-ft-deep, $50,000 to $500,000. 

9. April 25 0200 GMT, St. Louis, St. Louis County, Missouri: hail 
l-in.-diam (golf ball), thunderstorms. 
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10. April 25 0200 GMT, Cedarville, Crawford County, Arkansas: funnel 
cloud V « northeast, 

11. April 25 0549 GMT, Conway, Faulkner County, Arkansas: funnel 
cloud V - east. 

12. April 25 0500^-0555 GMT, Waldron, Scott County, Arkansas: thunder- 
storm, severe lightning. 

13. April 25 0608 GMT, Walnut Ridge, Lawrence, County, Arkansas: 
funnel cloud 1? = east. 

14. April 25 0630 GMT, Cotton Plant, Woodruff County, Arkansas; 
funnel cloud 7 - east. 

15. April 25 0637 GMT, Forrest City, St. Francis County, Arkansas: 
funnel cloud ? = east. 

16. April 25 0620 GMT, Craighead County, Arkansas: line of thunder- 

storms 7 = southeast, wind, $50,000 to $500,000. 

17. April 25 0645 GMT, Turrell, Crittenden County, Arkansas: tornado 

? = southeast, 40 yd/0.1 mi. 

18. April 25 0800-1200 GMT, southeastern Illinois: rain, earthen 

dam collapse. 

19. April 25 0200-0800 GMT, southern one-third of Illinois: hail, 
thunderstorm, $5,000 to $50,000. 

20. April 24-25, southern Indiana: heavy rain, $500,000 to $5,000,000. 

21. April 24-25, west and south central Kentucky: thunderstorms, 

rain, flooding. 

22. April 25 1030 GMT, Water Valley, Itawaba County, Mississippi: 
winds . 


23. April 24 1000 GMT, northern Alabama: thunderstorms ? ■ east, 
tornado, $50,000 to $500,000. 

24. April 24 1630 GMT, Weakley County, Tennessee: tornado, 

200 yd/0.5 mi . 

25. April 24 2000 GMT, Macon County, Tennessee: winds, $5,000 to 
$50,000. 

26. April 24 2230 GMT, Cumberland County, Tennessee: tornado, 

440 yd/1.5 mi., $50,000 to $500,000. 

27. April 25 0200 GMT, Tipton County, Tennessee: tornado, 

400 yd/6. mi,, $5,000 to $50,000, 

28. April 25 0720 GMT, Crockett County, Tennessee: tornado, Cairo, 
Nance, Quincy communities, 440 yd/1 1 mi., $50,000 to $500,000. 

29. April 25 0740 GMT, Humboldt and Gibson Counties, Tennessee: 
winds, $50,000 to $500,000. 

30. April 25 0800 GMT, Carroll County, Tennessee: winds, $5,000 to 


